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TECHNICAL NOTE 3080 

MEASUREMENT AND ANALYSIS OF WING AND TAIL BUFFETING LOADS 
ON A FIGHTER-TYPE AIRPLANE 
By Wilber B. Huston and T. H. Skopinski 

SUMMARY 


The buffeting loads measured on the wing and tail of a fighter- 
type airplane during 194 maneuvers are given in tabular form, along 
with the associated flight conditions. Measurements were made at 
altitudes of 30,000 to 10,000 feet and at speeds up to a Mach number 
of 0.8. Least-squares methods have been used for a preliminary analy- 
sis of the data. 

In the stall regime, the square root of the dynamic pressure was 
found to be a better measure of the load than was the first power. The 
loads measured in maneuvers of longer duration were, on the average, 
larger than those measured in maneuvers of short duration. Consider- 
able load alleviation was obtained by a gradual entry into the stall. 

In the shock regime, the magnitude of the load at a given speed and 
altitude was determined by the extent of the penetration beyond the 
buffet boundary. For a modification of the basic airplane in which 
the wing natural frequency in fundamental bending was reduced from 
11.7 to 9.3 cps by the addition of internal weights near the wing tip, 
a 15-percent decrease in wing loads and a similar percentage increase 
in tail loads resulted. 

The loads on a simplified wing buffeting model are examined on the 
assumption that buffeting is the linear response of an aerodynamically 
damped elastic system to an aerodynamic excitation which is a stationary 
random process. The agreement between the results of this analysis and 
the loads measured in stalls is sufficiently good to suggest the exami- 
nation of the buffeting of other airplanes on the same basis. 

INTRODUCTION 

An early investigation of buffeting which utilized the North 
American F-51D airplane (ref. l) provided basic information on the flight 
conditions under which buffeting was encountered and provided measure- 
ments of the magnitude of the buffeting loads on the horizontal tail. 
Speed and altitude were shown to be primary variables, and the load data 
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were reduced to dimensionless coefficient form by means of the product: 
Dynamic pressure x Tail area. It was hoped that such a buffeting- load' 
coefficient might be applicable to other airplanes, but the assumption 
that a form of coefficient common in steady-state aerodynamics would 
be applicable to a dynamic phenomenon was recognized as requiring 
further investigation. 

Since the completion of the tests of reference 1, a number of other 
experimental flight and wind-tunnel studies have been conducted. The 
effects of airfoil section and plan form on buffeting have been investi- 
gated. Buffet boundaries of a number of specific airplanes have been 
obtained. In several instances wing and tail loads have been measured 
during buffeting with special research airplanes. An analytical approach 
has also been made to the buffeting- loads problem, based on methods 
developed in the study of stationary random processes (see ref. 2). 

Upon completion of the tests of reference 1, plans were made to 
extend these tests of the same airplane to measure wing loads and tail 
loads simultaneously during buffeting, and, at the same time, to measure 
the effect of maneuver rate and the effect of penetration beyond the 
buffet boundary. In addition, the altitude coverage was to be improved 
in order to resolve more clearly the effect of this variable and, since 
it was thought that structural frequency might also be a significant 
variable, provision was made to modify the wings for several tests in 
order to measure some buffeting loads with a reduced wing frequency. 

The purpose of the present paper is to report the results of these 
extended flight tests and, especially, to report the magnitude of the 
buffeting loads measured. The basic load data involving 19A runs are 
given in tabular form, together with associated flight conditions. The 
results of preliminary studies which illustrate certain trends in the 
data are also given, but this analysis is not intended to be definitive. 
Although the present tests do not cover either the configurations or 
the speed range of greatest current interest, some of the variables are 
covered more extensively than in other tests. Stall buffeting, in 
particular, which will probably be common to all airplanes whatever the 
configuration, is extensively covered, and it is believed that all the 
data may be of value to those who are interested in the prediction of 
buffeting loads. The results of an analytical study in which the methods 
of generalized harmonic analysis are applied to a simplified wing buf- 
feting model are given in an appendix. 


SYMBOLS 


(Note: Symbols used only in appendixes are defined where they 

occur. ) 
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A 

A, B 
a, b 
b 

( C ki)eff 


C N 



k 

L 

AL 

M 

n 

P 

1 


aspect ratio, b/c 

constants used in tail-load equations 
constants used in wing- load equations 
wing span, ft 

effective slope of lift curve for damping of small oscilla- 
tions of a stalled wing in first bending mode 

airplane normal- force coefficient, nW/qS 

mean-square value of coefficient of section-normal- force 
fluctuations in buffeting ’ 

average wing chord, S/b 

frequency, cps 

pressure altitude, ft 

wing stiffness, lb/ft 

root structural shear load due to buffeting, lb 

amplitude of maximum root-structural-shear fluctuation due to 
buffeting encountered during run, lb 

Mach number 

normal load factor 

penetration beyond buffet boundary (defined in eq. (13) ) 
dynamic pressure, lb/sq ft 


r coefficient of linear correlation 

S area, sq ft 

s standard error 

t time, sec 

^load time between onset of buffeting and occurrence of measured 
load AL 
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V true airspeed, ft/sec 

W airplane weight, lb 

a angle of attack, radians 

(o circular frequency, 2 *f, radians/sec 

e residual, that is, a measured value minus a calculated value 

Subscripts: 

av average over class 

B onset of buffeting 

BB buffet boundary 

E end of buffeting 

L left 

max maximum 

n natural 

R right 

T tail 

W wing 

Mean values are designated by a bar (as c n ^); time differentiation 
by a dot (as a). 


AIRPLANE AND INSTRUMENTATION 
Airplane 


The airplane used for the present tests was the same North American 
F-51D airplane with heavily reinforced horizontal tail, fuselage, and 
wing used for the investigations reported in references 1 and 3 . A 
three-view diagram of the test airplane is shown in figure 1 ; a photo- 
graph, in figure 2 . 
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The airplane is equipped with a Packard V-165O-7, 12 - cylinder engine 
and a 4 -b laded Hamilton Standard Hydroraatic Propeller, 11 feet 2 inches 
in diameter. The propeller-to-engine gear ratio is 0 . 4-79 to 1 . Geo- 
metrical data for the airplane are listed in table I. The natural 
structural frequencies of various components as determined by ground 
vibration tests are listed in table II. In this table two sets of 
values of wing natural frequency are shown. One set applies to the 
basic airplane configuration and to the greater portion of the tests 
reported herein; the other set applies to the modified airplane, that 
is, the airplane with 100-pound weights added internally near the wing 
tips in order to lower the wing natural frequency in the fundamental 
bending mode from 11. 7 to 9-3 cps. 


Instrumentation 

Standard instrument s . - Impact pressure, pressure altitude, and 
normal acceleration were measured as functions of time with standard 
NACA recording instruments. The airspeed head was mounted on a boom 
extending 1.2 chords ahead of the leading edge of the wing near its 
right tip, and the NACA airspeed-altitude recorder was located near the 
boom to minimize lag effects which are believed to be negligible for 
the rates of change of altitude or airspeed encountered . The airspeed 
system was calibrated for position error up to a Mach number of 0.78; 
this calibration made possible the determination of the flight Mach 
number to within ±0.01. 

Airplane normal force was measured with an accelerometer mounted 
near the airplane center of gravity. The sensitive element had a 
natural frequency of 16 cps and was air damped. The damping was adjusted 
to 0.6 of critical at sea level, except during the tests with the modi- 
fied wing, when the damping was changed to 0.6 of critical at a pressure 
altitude of 30,000 feet. 

Strain-gage installation .- Measurements of structural shear on the 
wing and horizontal tail were made by means of wire resistance strain 
gages wired in four- active- arm bridges and attached near the roots of 
the principal structural members. Shear bridges were attached to the 
spar webs and bending- moment bridges, to the spar flanges. The entire 
installation was calibrated by established methods. (See ref. 4 .) For 
the shear on a wing panel, this calibration resulted in two combined 
strain-gage channels. One of these combined channels was principally 
sensitive to shear and secondarily sensitive to bending moment; the 
other channel was primarily a measure of bending moment, and secondarily 
sensitive to shear. The outputs of these two channels, recorded as a 
function of time on a multiple-channel recording oscillograph, could be 
combined numerically to obtain the wing-panel structural shear. The 
shear on the left and right panels of the horizontal stabilizer was 
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obtained from the outputs of the left and right combined strain-gage 
channels which were sensitive to shear. This strain- gage system repre- 
sents an improvement over that used in reference 1. 

The recording oscillographs used employed galvanometer elements 
with a natural frequency of 100 cps which were damped to about 0.6 of 
critical damping. This combination of damping and natural frequency 
insured an approximately linear response for the buffeting frequencies 
expected. Special care was taken to balance the galvanometer elements 
so as to keep any possible acceleration effects within the reading 
accuracy. Variations in sensitivity due to voltage changes were elimi- 
nated by provision of a calibrate signal on the record for each run, 
and the stability of the strain-gage installation was checked at 
intervals by application of known loads to the wing and tail. The 
overall experimental error in incremental values of wing root shear 
obtained from the strain-gage— oscillograph system is estimated from 
the calibration as less than ±130 pounds; while for the incremental 
values of shear on the right and left horizontal stabilizer the esti- 
mated error is of the order of ±80 pounds. 


TESTS 


All tests were made with the airplane in the clean configuration, 
and the power setting, at low Mach numbers, was that required to attain 
level flight at the altitude of test. In tests at Mach numbers greater 
than the level-flight capabilities of the airplane, normal rated power 
was used. Of a total of 194 runs in which buffeting was measured, 150 
were made with the basic airplane and 44 with the modified airplane. 

With the basic airplane, gradual turns to the stall were performed 
at nominal test altitudes of 30,000, 25,000, 20,000, 15,000, and 
10,000 feet. Pull-ups were performed at 30,000, 25,000, and 20,000 feet. 
The range of Mach numbers covered was 0.34 to O.792 at 30,000 feet and 
0.23 to 0.4l at 10,000 feet. 

With the modified airplane, the added wing-tip weights introduced 
local stress concentrations which restricted the maximum allowable load 
factor for buffeting flight to 4, and limited the maneuvers to pull-ups. 
With the airplane at 30,000 feet, buffeting cannot be obtained at speeds 
between M = 0.54 and M = 0-73> without exceeding the limit load fac- 
tor of 4; while at 10,000 feet, buffeting is not encountered at speeds 
between M = 0.32 and the maximum permissible diving speed which for 
the standard North American F-51D airplane is a true airspeed of 537 mph. 
For the modified airplane, buffeting was, therefore, obtained by per- 
forming pull-up maneuvers at 30,000 feet and 10,000 feet at speeds 
limited by the foregoing considerations. 
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METHOD OF OBTAINING DATA 


The procedure and definitions used in presenting the results of this 
investigation are best illustrated by referring to the typical time- 
history records shown in figure 3- The accelerometer record (fig. 3(a)) 
was used to establish the time for the beginning tg and end tjj of 

buffeting, as well as the duration of buffeting. These values were 
obtained simply by observing the point at which there was a distinct 
change in the character of the accelerometer trace . The airplane normal- 
force coefficient Cn was obtained from the accelerometer and airspeed 
records . Values of Cn during buffeting were based on a mean line 
faired through the fluctuations of the accelerometer record. The air- 
plane normal-force coefficients at the beginning and end C^g of 

buffeting were determined and corresponding values of Mach n umb ers 1% 
and Mg were also noted. In determining all values of airplane normal- 
force coefficient, the value of airplane weight W used for each run 
was the take-off weight corrected for the fuel consumed prior to the 
start of the run. The maximum rate of change of airplane normal load 
factor n prior to the onset of buffeting was determined for each run, 
as in figure 3(a), and the maximum rate of change of angle of attack 
per chord traveled ac'/v was estimated from h on the assumption that 
the speed remains constant and 


. dC]\j/ dt nW 

OC = — 

dC^/da qS^dCjj/da) 

and hence that 


dc ^ n C ^B c 
V dCjj/da V 


In this relation, a nominal value of 5*3 was used for dC^/da. 

A typical .oscillograph record for obtaining wing and tail loads is 
shown in figure 3(b). The six traces identified with numbers in this 
figure were employed. Traces 1 and 2 are measures of root shear on the 
right and left horizontal tail, respectively. Root shear on the left 
wing panel is measured by a combination of the deflections of traces 15 
and 17; on the right wing, by a combination of traces 5 and l6. Buf- 
feting loads, which are incremental loads, were determined from the peak- 
to-peak deflections of these traces (designated S]_, etc., in fig. 3(t>)). 
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The buffet-load values AL reported for a run are one-half of the 
largest peak-to-peak fluctuation in each of the four loads encountered 
during that run. The time of each load maximum was recorded and is 
reported as the incremental time Ati oa a following the onset of buf- 
feting. Through use of a timer common to the standard flight instru- 
ments, values of M, Cjj, and q corresponding to each buffeting load 

were determined. 

RESULTS 

Buffet Boundary 


The data acquired in the present investigation of the basic air- 
plane are incorporated in table III. For the modified airplane the data 
are included in table IV. Tables Ill(a) and IV(a) deal with the oper- 
ating conditions under which buffeting was first encountered and under 
which it ended. In addition to the numerical data, a pilot's note 
column is included. In most instances the pilot estimated the intensity 
of buffeting, in one of four categories: very light, light, moderate, 

or heavy. These comments have been designated by the letters vZ, Z, 
m, and h. The pilot's notes on the direction of the roll-off after 
the stall are also included, left and right roll being desi gnat ed by L 
and R, respectively, while no roll is indicated by N. 

The flight conditions for the onset and end of buffeting given in 
tables Ill(a) and IV(a) are summarized in plots of airplane normal-force 
coefficient against Mach number in figures 4 and 5 , respectively. In 
figure 4(a) a buffet boundary for the onset of buffeting is also shown 
and two labels Stall regime and Shock regime are included. These 
labels denote speed regimes in which the flight characteristics of the 
airplane differ, and thus speed regimes in which the buffet boundary 
was obtained in different ways. For Mach numbers below about O. 65 , 
buffeting was usually encountered in an accelerated stall maneuver; 
a maximum value of airplane normal-force coefficient was reached; and 
controlled flight at still higher load factors was not then possible. 

In this stall regime the value of Cft-g for the onset of buffeting 

varied with Mach number and also was generally higher in pull-ups than 
in turns. The increase can be associated with the abruptness of the 
stall entry, as measured by the largest value of dc/V reached prior 
to the onset of buffeting. The buffet boundary shown for the stall 
regime in figure 4(a) was obtained from faired cross plots of C$g, M, 

and dc/V, greatest weight being given to the data for 30,000 feet, and 
corresponds at each Mach number to the value of Cjfg for dc/v = 0. 

The difference between this boundary and the actual Cjjg at' the onset 
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of buffeting is plotted as a function of oc/v in figure 6 for the data 
from altitudes of 30,000, 20,000, and 10,000 feet. The increment in 
normal-force coefficient is analogous to the increment in the dynamic 
value of the maximum lift coefficient as compared to the static value, 
but, because of the approximate nature of the relation between acceler- 
ometer reading and rate of change of angle of attack, a more detailed 
study which might include the effects of Reynolds number has not been 
attempted. For this reason also, no attempt has been made to specify 
a variation of buffet boundary with altitude, although the possibility 
of such a variation is suggested by a comparison of the plots for 
30,000 feet and 10,000 feet in figure 6. 

For Mach numbers above about O. 65 , buffeting was encountered during 
diving turns or in pull-outs from dives. The onset of buffeting occurred 
at values of Cjj well below maximum lift, but controlled flight at 
normal-force coefficients well above the value for the onset of buffeting 
was feasible. The buffet boundary shown in figure 4(a) above M = 0.64 
was obtained by fairing through the observed values of Cjjg, greatest 
weight being given to the data for 30,000 feet. 

The buffet boundary of figure 4(a), based on data for the onset of 
buffeting, appears to define a transition from steady to unsteady 
phenomena. This boundary, which has been placed in figure 4(b) for 
comparison, does not appear to define the transition from unsteady 
back to steady conditions. The data for the end of buffeting represent, 
however, the flight conditions on final subsidence of oscillations in 
the structure. In the shock regime, when buffeting persisted to values 
of Cn below the buffet boundary and the return to level flight from 
the maximum load factor was rapid, the persistent fluctuations appeared 
to differ in character from the rest of the record, resembling the 
subsidence of a damped oscillation from which the excitation has been 
removed. When the approach to the boundary was at a slow rate (generally 
accomplished by a loss of speed at nearly constant load factor), the 
end of buffeting occurred as the boundary was crossed. The buffet 
boundary above M = O .65 as defined by the onset of buffeting may, 
therefore, represent a distinct boundary below which a buffeting excita- 
tion is not present. 

In the stall regime, values of in almost all instances are 

below the buffet boundary. Although the persistence of structural 
oscillations may be a factor in this case also, the character of the 
fluctuations indicates that buffeting, once encountered, is maintained 
to values of Cjj reached in the stall recovery which are well below 
the buffet boundary. 

The buffet boundary for the basic airplane, figure 4(a), has been 
plotted in figure 5 ( a ) for comparison with the data for the modified 
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airplane. The boundary for the basic airplane appears to represent the 
modified airplane reasonably well. The two points for C^g at the 

lowest Mach numbers are for maneuvers at 10,000 feet, and may represent 
a Reynolds number effect, but enough data to establish a consistent 
trend are not available. 


Wing and Tail Buffeting Loads 

The wing buffet loads associated with the runs of table Ill(a) and 
IV(a) are given in tables Ill(b) and IV (b); the tail buffet loads are 
given in tables IIl(c) and IV(c). There is also listed a quantity AC^, 
the penetration beyond the buffet boundary in terms of mean airplane 
normal-force coefficient, used in the analysis of sbme of these data. 

The wing and tail buffet -load values for the basic airplane given 
in tables Ill(b) and IIl(c) are shown in summary form in figures 7 and 
8 ; the data for the modified airplane are shown in figures 9 and 10 . 

In these figures the variation of the loads on the left and right sur- 
faces with Mach number is shown for each of the nominal test altitudes. 
Turns are distinguished from pull-ups . 

In the absence of any accepted theory relating the magnitude of the 
loads in buffeting to the flight conditions and the characteristics of 
the structure, the analysis of the load data of tables III and IV has 
necessarily been of a somewhat qualitative nature, involving both 
general regression studies and the fitting of regression equations to 
the data by means of least-squares methods. The results of this study 
are incorporated in the following section. 


ANALYSIS AND DISCUSSION OF BUFFETING- LOAD DATA 


When the buffeting-load data of tables III and IV are plotted 
against Mach number for different altitudes, the large amo unt of scatter 
in, for example, figures 7 and 8 makes it difficult to assess the effects 
of both speed and altitude and suggests that other factors may be signif- 
icant. As shown by the difference between the data for turns and pull- 
ups in figure 7 (a), one such factor is the abruptness with which the 
stall is entered. A number of studies have been undertaken in attempts 
to identify other significant parameters. In these studies use has been 
made of the usual methods of regression analysis, including correlation 
studies, graphical studies, and the fitting of regression equations by 
least-squares methods. The form of these equations was inferred from 
the graphical studies or in some instances could be based on analytical 
results. In these studies the loads measured in stalls were found to 
follow a somewhat different pattern from those measured in the shock 
regime . 
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As a preliminary to analysis of the load data, a considerable 
simplification was effected on the basis of plots of left wing load 
against right wing load and left tail load against right tail load 
shown in figures 11 and 12, respectively. The coefficient of correla- 
tion shown in these plots, of the order of r = 0.9, can be regarded 
as a measure of common causes and suggests that the factors which pro- 
duce loads of a given size are in general common to the left and right 
wing panels, or left and right tail surfaces. On this basis, the mean 
value Aly of the two wing-panel loads measured in a run was taken as 
representative of the wing loads encountered during that run; that is, 
the mean wing load Aly = 0.5(AI^k + Alyp) and a similar mean tail 
load AI/p = 0.5(ALpl + AI/jr) were used to represent the loads in each 
run. 


A scatter diagram of AI# against AL? is shown in figure 13 . 

The value of the coefficient of correlation, 0.7, suggests a larger 
degree of independence between wing and tail loads than is the case 
for the left and right wing or tail surface. On this account, analysis 
of the wing and tail loads was carried out independently. 


Regression Analysis 

When dealing with quantities of data, the interrelation of more 
than two parameters cannot ordinarily be shown in a simple plot, but 
the effect of a given independent variable can be investigated if the 
data are grouped by classes of this variable and the average values of 
the dependent variable (in the present case the load AL ) are computed 
for each class. Provided that each class constitutes a similar sample, 
the effect of other independent variables on the load may thus be sup- 
pressed, or averaged out, and the variation with the independent variable 
of interest established. The grouping and averaging may then be repeated 
for other variables. Such an analysis is, of course, somewhat qualita- 
tive, and it may be difficult to show the effect of a secondary variable 
in the presence of a large primary effect. 

In the study of loads measured on the basic airplane, the variables 
investigated for runs in which the stall was reached include dynamic 
pressure q and the length of time spent in buffeting At. Also inves- 
tigated was the effect of the abruptness of the stall entry. For this 
investigation the value of ac/v was used as a measure of the abruptness 
of the entry in both turns and pull-ups. For buffeting encountered in 
the shock regime, the variables investigated include the dynamic pressure 
and the increment in normal-force coefficient beyond the buffet boundary 
at which the load AL was measured. The trends shown by this study 
for both the stall regime and the shock regime are presented in the four 
parts of figure 14. 
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Load, trends in stall regime .— Stall buffeting in the present study 
occurs at Mach numbers below a value estimated, as O.65 ± 0.01. All 
runs in table Ill(a), therefore, for which % < 0.64 and for which 
values of ac/v and At could be established were included in the 
stall analysis. For each of the 91 runs thus available, the wing- load 
value Alty and the tail- load value ALp were used, together with the 
mean of the dynamic-pressure values, tables Ill(b) and IIl(c). 

The average variation of wing load with q is shown in figure 14(a). 
For this plot, the values of ALy were grouped into eight classes, 
according to the value of q; the plotted variable (AL^) av is the 

average of the loads Aly in each class. For the stall regime, the 
dynamic pressure increases by roughly a factor of 4 (i.e., 42 to 
180 lb/sq ft) while the average load increases by a factor of only 2 
(i.e., 500 to 1,000 pounds), an increase which is roughly proportional 
to the square root of q. The dynamic pressure is thus revealed as a 
major parameter in stalls, but the relation to load appears to be 

Alty oc ^q rather than Al^ oc q. This proportionality is used to examine 
the variation of wing loads in stalls with maneuver abruptness and with 
time spent in buffeting in figures 14(c) and 14(d), respectively, where 

plots of (aI^ / / q) av against ac/V and At are shown. An alleviating 
effect on load associated with a gradual stall entry is indicated since, 
at ac/v « 0, the loads (expressed as AL^ /i/ q) are as much as 40 percent 
less than the loads measured in more abrupt maneuvers where dc/V « 0.008 
radian per chord. The alleviation is indicated in figure 14(c) to be 
somewhat exponential in character. With regard to time spent in buf- 
feting, figure 14(d) suggests that on the average the maximum load 
encountered during buffeting increases with the total duration of 
time At spent in buffeting. From periods of less than 1 second to 
periods of 4 to 5 seconds , the increase is of the order of 90 percent 
but does not appear to be linear. 

The trends shown qualitatively in figures 14(a), 14(c), and 14(d) 
suggest a number of equations which can be written relating wing load 
to various combinations of the variables ■ representing speed, altitude, 
time, maneuver abruptness, and structural frequencies. Among the equa- 
tions investigated for the wing loads in stalls were the following: 


ALtf = 

(1) 

Altf = a 2 q 

(2) 


» 


Alty = 


(3) 
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ALtf = a^q log e (f n At) 

AI* - (a, ♦ v e^/b.OOM^ 


(4) 

(5) 


at ( . -u -a-c/O.OOMA / — 7- s 

% = V a 6 + b 6 e y/q. logg ( f n At) 


( 6 ) 


The values of the arbitrary constants in equations (l) to (6) can be 
obtained by fitting the equations to the experimental data. An advantage 
of the least-squares method of fitting lies in the ready availability of 
precision measures for the constants and of the standard error of esti- 
mate of the equation. (For convenient reference, definitions of terms 
and a summary of least-squares procedures as used in the present inves- 
tigation are included in appendix A. ) The results of the least-squares 
analysis of the wing loads in stalls are given in table V which shows 
the equations, the sums of the squares of the residuals, and the 
standard errors of estimate of the equations, together with the numeri- 
cal values of the constants and their standard errors of estimate. 

Equation (l) is of chief interest for comparison purposes. The 
value aj_ = 7^9 pounds in table V is the mean of the 91 values of AL^ 

being analyzed. The standard error of estimate, 255 pounds, is in a 
sense a measure of the error involved in the simple assumption that the 
data on the wing buffeting loads in stalls can be represented by this 
mean value . 

Equations (2) and (3) represent the combined effect of speed and 
altitude. Equation (2) is analogous to the dimensionless coeffi- 
cient Cg = ^ which parallels the usual coefficients for steady aero- 
dynamic forces and which has been much used in buffeting studies . 

Equation (3), which was proposed in reference 5 and which also follows 
from the analysis in appendix B, represents the combined effect of an 
aerodynamic excitation and an aerodynamic damping. The standard errors 
of estimate for these equations, 293 pounds and 226 pounds, appear to 
indicate that q is not so good an indicator of the size of the load 
as is the mean value, while y/q is better than the mean. A dependency 
of load on the square root of the dynamic pressure is also in line with 
the indications of figure 14(a), for stalls. Superiority of the square 
root of the dynamic pressure (as a measure of buffeting) as compared to 
the first power indicates that in stalls at a given altitude the loads 
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would be directly proportional to the Mach number or the true airspeed, 
while at a given Mach number (or airspeed) the loads would vary directly 
as the square root of the atmospheric pressure (or density). The linear 
trend with Mach number revealed by the least-squares analysis is recog- 
nizable in the data of figure 7 for stalls when, as for example in fig- 
ure 7(a ) > enough runs are available to give a representative distribution 
of the time spent in buffeting and the abruptness of the stall entry. 

The trend with pressure at a given Mach number is less evident, but for 
a pressure change from 628 lb/sq ft at 30,000 feet to 1,455 lb/sq ft at 
10,000 feet, the corresponding load increase is clearly less than the 
ratio of the pressures (2.32) and more nearly the square root of the 
pressure ratio (1.52). 


With regard to equation (4) in table V, it would ordinarily be 
expected that for a process in which random factors play a part, the 
probability of occurrence of a given value is higher for a large sam- 
ple than for a small one. The indication in figure 14(d) that larger 
loads are encountered in stalls of longer duration is qualitative con- 
firmation of this expectation. For a stationary random process, as 
outlined in appendix B, analytical results are available for determining 
the probability that a given peak value will occur once in a time At. 
These results lead to equation (4), and the standard error of estimate, 
206 pounds, represents an improvement over equation (3). In determining 
the value of a^, the value of the frequency of wing fundamental bending 
(ll.7 c P s j table II) was used for fn« This frequency is the one most 
often observed in the wing-shear strain-gage records. 


The roughly exponential trend of the variation of 




with 


'av 

oc/V indicated in figure 14(c) suggested the form be“® <? / V X Constant 
as a measure of the effect of maneuver abruptness on the loads in stalls. 
This form is purely empirical and was adopted simply to account in an 
approximate way for the observed trend in the data. Although a value 
of the exponential constant could have been determined by nonlinear 
regression methods, reference 6, the iterations required make the deter- 
mination much more laborious than the evaluation of the constants of 
the linear variations. Preliminary investigations having indicated a 
value of approximately 0.004 for the constant; this value was used in 
equations (5) and (6). In comparing equation (5) with equation (3) or 
equation (6) with equation (4), the relative magnitudes of the standard 
errors of estimate indicate a significant improvement resulting from 
inclusion of a measure of the maneuver abruptness. The relative values, 
of ag and bg (that is, 65.6 and -31.6) indicate that a load allevia- 
tion of about 50 percent could be obtained by a gradual stall entry. 
Although the physical basis for this alleviation is not understood, it 
may be associated with a less completely developed stall in the slower 
maneuvers resulting from a less abrupt flow breakdown. A brief study 
of the correlation between the duration and abruptness of the maneuvers 
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included in the analysis indicates that the larger loads in abrupt 
maneuvers were not explainable on the basis of stalls of longer dura- 
tion, but the magnitude of the effect of abruptness indicates that this 
factor warrants further examination and should not be ignored in other 
studies of wing buffeting loads in stalls. 

For the analysis of the tail loads in stalls, the equations 
examined include the following: 


= 0 


ATrp = Ay 

(7) 

ALp = Agq 

( 8 ) 

1 

II 

(9) 

AI? = Aio/q log e (f n At) 

( 10 ) 

A _ -6£/0.00kv\ r- 

- (A 1± + B i:L e J j/q 

( 11 ) 

-oc/O.OOl+v'N / 7 : 

5 + B 12 e Jyq. log e (f n At) 

( 12 ) 


The results of the least -squares analysis shown in table VI are 
for the same 91 maneuvers used in the wing- loads study. The form of 
equations ( 7 ) to (12) parallels the form of the equations used in the 
wing- loads study. Because of the empirical nature of the abruptness 

alleviation expressed by the term e -<^/O.OdkV , the wing chord and the 
constant 0.004 were retained in the tail- load calculations. The wing 
natural frequency was also retained in the expression log e (f n At). 

Comparison of the standard errors of estimate of the equations of 
table VI indicates the pertinence of the square root of the dynamic 
pressure, the duration of the stall, and the abruptness of the maneuver. 
The load alleviation obtainable by a gradual stall entry appears to be 
even greater than in the case of the wing loads . 
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Load trends in shock regime .- Buffeting at the Mach numbers of the 
shock regime was, for the present airplane, encountered under transient 
conditions in diving turns and pull-ups. In some instances so much 
speed was lost during a maneuver that buffeting originally encountered 
at a Mach number of 0-7 ended at Mach numbers of 0.62 or O.63 with a 
typical stall recovery. In order to assure a homogeneous class of data, 
the 26 runs selected as representative of the shock regime were those 
in which the maximum buffeting load was encountered at Mach n umb ers 
above 0.68, as shown by the Mach numbers of tables Ill(b) and IIl(c). 

A plot of values of (ALw)av against q for these maneuvers, fig- 
ure 14(a), appears to indicate a different trend with dynamic pressure 
in the shock regime than in the stall regime. One reason for the 
apparent trend with q is found in an examination of the variation of 
load with penetration beyond the buffet boundary. At a given Mach num- 
ber, increasing penetration beyond the' buffet boundary results in 
increased amplitude of load fluctuation, but the rate of increase of 
load with penetration varies with Mach number. These trends for the 
wing loads in the shock regime are illustrated in figure 15. 

Figure 15(a) shows the wing- load values AL^ plotted on a diagram 
of the variation of with Mach number. In each symbol is a numeral, 

indicating the value of AI*/ in hundreds of pounds. Also shown is the 
buffet boundary for the shock regime from figure 4. In general, smaller 
loads occur near the buffet boundary, larger loads at values of 
farther removed from the boundary. Figure 15(b) is a plot of load 
against the difference ACjj = Cjj - for Mach numbers of approxi- 

mately 0.7 and 0-75* The linear dependence of load on ACjy is evident, 
but the slope dALw/dAC^ decreases as M increases. 

Shown also in figure 15(a) is a line marked This curve of 

maximum normal-force coefficient was estimated from a study of recent 
wind-tunnel data on C Nmay since specific data for the North American 
F-5ID are not available. If the penetration beyond the buffet boundary' 
at each Mach number is expressed as a ratio denoted by P where 


P = 


C N - %b 

^max 


(13) 


the Mach number dependence of the slopes in figure 15(b) is accounted 
for. A plot of (Alw)av against P is shown in figure 14(b). The 
variation of (Alty)av with p appears to be linear for the range of 
flight-test data available; the strong dependence on P effectively 
masks any dependence on q in figure 14(a). 
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The equations investigated for wing loads in the shock regime were 


My = a l4 

(14) 

My = a 15 q 

(15) 

% = a i6^ 

(16) 

ALtf — a ±r jP 

(IT) 

AI^ = a lQ Pq 

(18) 

£*w = a i9 P /T 

(19) 


The results of the least-squares analysis are given in table VII. 

For the tail loads in the shock regime the equations investigated 
were similar to those for the wing loads, that is. 


AD? = A 2 o 

(20) 

ALp = 

(21) 

ADp = A22 \/q 

(22) 

ALp = A2^P 

(23) 

ALp = A24Pq. 

(24) 

ALp = A2^P\/q 

(25) 


The results of the least-squares treatment are shown in table VIII. 
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For both wing loads and tail loads in the shock regime, the values 
of the standard errors of estimate show that neither q nor i/q is as 
good a measure of the load as the average value, although /q is some- 
what better than q. Inclusion of the penetration in the analysis 
through the parameter P (eqs. (17), (l8), (19), (23), ( 2 k), and (25)) 
results in values of the standard error of estimate which are clearly 
very much lower than the values for the means (eqs. (l4) and (20)). 
Between equations involving P, Pq, and P^j, the indications are not so 

clear. For wings equation (19), AL^ = a-^P/q, has the smallest stand- 
ard error of estimate, while for tail loads equation (23), ALq = A23P, 
has the smallest standard error of estimate. The lack of a clear indi- 
cation about the effect of q in the shock regime may be in part the 
result of the relatively small number of points and the limited range 
of altitudes that are available at a given Mach number. A further con- 
tributing factor may lie in the random character of the buffeting proc- 
ess as discussed in appendix B. The strong dependence of resultant 
loads on penetration, coupled with the transient character of the maneu- 
vers at speeds above the maximum speed in level flight, would require a 
more detailed analysis including perhaps not only the extent of penetra- 
tion but also the length of time spent at or near any given value of 
penetration. Since the standard errors of estimate for equations (23), 
(24), and (25) are so nearly the same, it will be assumed that the 
variable is also applicable to the tail loads in the shock regime. 


Load Equations of Best Fit 

Wing loads .- The summary of the regression analysis of the wing 
loads measured in the present tests, tables V and VII, indicates that 
the best fit is obtained with equations (6) and (19). These equations 
may be written in terms of the values of the regression coefficients 
as, for the stall regime, 


Aiy = 65.6 ± 3.8 - (31.6 ± 5 .io e -<Wo.oo4v ^ loge ( 11>7 At ) 


(26) 


and, for the shock regime, 


Altf = (153.5 ± 6 .4 )P j/q 


(27) 


In figure 16 a comparison is made of the variations of wing load given 
by equations (26) and (27) with the effects of q, maneuver abruptness, 
stall duration, and penetration shown in figure 14. The data points of 
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figure 16 are reproduced from figure lit-. Shown in each part of the fig- 
ure are the mean values of the "suppressed" independent variables. For 

the stall regime, these values (ac/v) = 0.00193 radian per chord and 
At = I.78 seconds have been substituted into equation (2 6) in order to 
show in turn the variation of (AIy) av with q, figure 16(a), the 

variation of (AL^//q^ av with ac/V, figure l6(c), and the variation 

of (aI w/^0 av ^th At, figure 16 (d). In the shock regime, the average 

value of q has been substituted into equation (27) to show the trend 
of (AL^) av with penetration P. (See fig. 16 (b).) Since the trend 

of load with q in the shock regime has been obscured by the large 
range of values of penetration P, no comparison is shown in figure 16 (a). 
The agreement between the points representing average trends and the 
dependency on \/q and At in equation (26) is substantial and suggests 
the validity, at least for the present airplane, of the physical con- 
cepts represented in the form |/q log e (f n At). The exponential character 
of the alleviation in load obtainable by a gradual stall entry, even 
though empirical, appears also to represent the trend in the experi- 
mental data. Since the effects of duration and abruptness can both be 
of the order of ±25 percent of the load for an average condition, the 
advisability of examining the buffeting of other airplanes on the same 
basis is indicated. 

The expression of the penetration beyond the buffet boundary by 
means of the ratio (C N - % B )/( C IW " Cn Bb) as in equation (13) is 

purely empirical, but over the range of flight-test data available 
appears to give a reasonably good fit to the data (fig. 16 (b)). The 
linear dependency of load on P assumed in the regression analysis is 
also empirical, and verification for large penetrations at Mach numbers 
above 0.70 is not feasible with the present airplane because of opera- 
tional limits. In particular, it is not known whether the loads for a 
stall at transonic speeds would be given correctly or whether, as at 
lower speeds, the abruptness of stall approach would be important; 
investigation with an airplane with wider operational limits is desirable. 

/ 

A comparison of the loads calculated by use of equations (26) and 
(27) with the measured loads on which the numerical values of the regres- 
sion coefficients are based is shown in figure 17. In each part of fig- 
ure 17 the line of exact agreement is the solid line with unit slope. 

The horizontal or vertical distance from any point to this line is the 
difference between the measured and the calculated load. Parallel to 
each line of exact agreement are two dashed lines, displaced by the 
amount of the standard error of estimate. In general, 68 percent of 
the measured values will vary from the calculated values by less than 
the amount of the standard error of estimate. The wing loads calculated 
from equations (26) and (27) when compared with the measured values 
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(figs. 17(a) and 17(b)) show generally good agreement. The measured 
wing loads are estimated to be in error by less than ±130 pounds, as 
compared with a standard error of estimate for equation (26) of 178 pounds 
and for equation (27) of 228 pounds. The fact that in the stall regime 
these two precision measures have roughly the same order of magn itude 
suggests that, with the present data, regression analysis can probably 
accomplish little more; in the shock regime, the larger standard error 
of estimate for equation (27) as compared to the error limits of the 
experimental data -may be a further indication of the need for a more 
detailed study than has been possible with the present data. 

Tail loads .- The summary of the regression analysis of tail loads 
measured in the present tests indicates that the best fit of the stall 
data (table VI ) is obtained with the equation 

Al/f = \kk.l ± 2.9 - (29.2 ± 4 .l)e" aC /°* 004 ^]/q log e (ll .7 At) (28) 

while the equation which is taken as representing the shock- regime data 
(table VIII) is 


ALj, = (75.2 ± 4 . 6 )P|/q (2 9 ) 

Loads calculated from these equations are compared in figures 17(c) and 
17(d) with the measured loads from which the regression coefficients 
were obtained. Since equations of the same form as the wing- load equa- 
tions give such a good fit, the possibility is indicated that the wing 
is a primary agency in determining tail loads. Since the response of 
the tail is primarily at a frequency corresponding to that of the fuse- 
lage in torsion, the wing may excite the tail through the fuselage. On 
the other hand, the standard errors of estimate for equation (28), 

135 pounds, and for equation (29), 174 pounds, are somewhat larger than 
the estimated experimental error (±80 pounds) and this difference, 
coupled with the correlation coefficient of 0.7 between tail and wing 
loads, indicates that one or more additional parameters may exist which 
are important in determining tail loads but which are not disclosed by 
the present investigation. The propeller slipstream may provide one 
such agency and the wing wake another, but, since instrumentation suit- 
able for the evaluation of such effects was not incorporated, the rela- 
tive contributions of the fuselage, the wing wake, and the propeller 
slipstream cannot be established. 
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Extension of Results 

Comparison of loads measured on basic and modified airplane .- The 
large amount of scatter in plots of buffeting load against Mach number 
in figures 7, 8, 9, and 10 makes difficult any simple determination of 
the effect of the added wing-tip weights on the magnitude of the buf- 
feting loads. Comparison of figures 7(a) and 9(a), for example, is 
inconclusive. The equations obtained in the analysis of the buffeting 
loads on the basic airplane have, therefore, been employed to extend 
the results obtained on the basic airplane to the analysis of the data 
for the modified airplane. For the stall regime, equations (2 6) and 
(28) have been used, modified only to the extent required to allow for 
the slightly reduced probability of encountering a given load in a 
given time since the wing frequency has been reduced. The equations are 


AIw = ^65.6 - 31 .6e'^°' 00 ^ V )^q loge(9-3 At) 


(30) 


ADf = (kk.l - 29. 2e aC /°’° (AV )/q loge(9.3 At) 


(31) 


In the shock regime equations (27) and (29) were used. Values of oc/v 
and At from table IV(a) were used with average values of q and ACfl 
from tables IV(b) and IV(c) to calculate values of AT^ and ALj. 


These calculated values are compared with the values measured in flight 
in figure 18 in which the solid lines are lines of exact agreement. 

As a measure of the effect of the reduced frequency on load, the average 
(AL)modified 

has been determined, by computing the value of k 


ratio 


(AJ j )basi C 
in the equation 


(AL)modified- k(AL)-b as i C 


(32) 


The values of k for the wing and tail in the stall regime and shock 
regime together with their standard errors of estimate are 


^wing, stall ” 0.90 t 0.03 

^wing, shock = O.7I ± 0.07 

k tail, stall = 1-25 ± 0.04 
^tail, shock = 1*10 ± 0.10 
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The dashed straight lines represented by these values of k are shown 
in figure 18. 

For the wing in the stall regime, the value of k indicates an 
average reduction of 10 ± 3 percent over and above the average reduc- 
tion of about k percent that would be expected because of the reduced 
probability associated with the frequency reduction. The est ima te of 
a 29 * 7 percent load reduction in the shock regime is somewhat less 
reliable than the 10-percent estimate since a smaller number of points 
is involved, but an overall reduction of something like 15 percent is 
indicated for the modified airplane. 

Comparison of the tail loads measured on the modified airplane 
with the loads calculated from the least-squares equations as shown 
in figures 18(c) and 18(d) indicates that the wing modification has 
increased the tail loads about 15 percent. In buffeting, the motion 
of the tail is primarily in an antisyrametrical mode at the natural 
frequency of the tail assembly as restrained in torsion by the fuselage, 
9.8 cps in table II. Since the addition of the wing-tip weights reduced 
the frequency of the wing in fundamental bending from 11. 7 to 9.3 cps, 
table II, wing buffeting of the modified airplane occurs at a frequency 
only about 0.5 cps removed from the tail buffeting frequency, whereas 
with the basic airplane, the difference is nearly 2 cps. The amplitude 
response of a simple system would be expected to be larger, the nearer 
the frequency of the excitation to resonance, and it is possible that a 
coupling exists between wing and tail vibration modes such that this 
simple explanation would be sufficient to account for the experimental 
results. If so, the importance of the fuselage as a coupling agent in 
the tail- load problem is indicated. 

Measured loads compared with results for simplified wing buffeting 
model . - In appendix B, an equation is developed which gives the form of 
the relation between pertinent structural and aerodynamic parameters 
and the mean-square value of the root-structural-shear fluctuations of a 
stalled wing under the assumption that such buffeting can be treated as 
the response of a damped linear elastic system to an aerodynamic excita- 
tion which is a stationary random process. The buffeting model considered 
is a simplified wing with one degree of freedom, fundamental bending, 
and the development parallels, in some respects, the study of the loads 
on a tail in a fluctuating airstream in reference 2. The development 
is tentative, since the assumption that stall buffeting is a normally 
distributed stationary random process has yet to be verified, but a 
comparison of the loads measured in the present study with the tenta- 
tive relation is of interest . 

A primary aerodynamic factor determining the magnitude of the buf- 
feting loads is the power spectrum of the aerodynamic excitation, 
denoted by the spectrum of the coefficient of sect ion-normal- force 
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fluctuations c n 2 (cu) in appendix B. Provided that this spectrum 
possesses certain general dimensional and frequency characteristics 
(especially a fairly constant level over a band of low frequencies), 
the details of the shape of the spectrum are of minor concern, but 
the mean-square value of the excitation c n 2 is of great importance. 

In appendix B, the scale factor in the power spectrum of the excitation 
is assumed to be the chord, the damping is assumed to be positive and 
aerodynamic, and the resultant equation for the root -me an- square shear 
at the root of a wing panel due to buffeting (eq. (B27)) is 




-. 1/2 



(33) 


In this equation the operating conditions of speed and altitude are 
included in the term /q; the geometry of the wing and its stiffness 
are included in the term in parentheses; while the ex cit ation and 


the aerodynamic damping are represented by the term c n 2 /(^i tL ^ ( 


eff* 


Little information is available about any spectrum of section normal 
force, or about the term which is an effective slope of the 

lift curve applicable to the aerodynamic damping of small bending 
oscillations of a stalled wing. Unpublished tests in the Langley 
2- by 4-foot flutter research tunne l on a stalled, rigid NACA 65AOIO 

airfoil have given values of / c ^ O.O7 over a range of angles of 


attack beyond the stall. Vibration tests of a similar stalled wing 
have indicated that over a wide range of reduced frequencies and angles 
of attack the aerodynamic damping is of the same order of magnitude as 
that indicated by the two-dimensional slope of the lift curve, that is, 
(Cl^eff 5,5 2n . Using these two results as a guide to order of magnitude 
gives a value 


c n‘ 


W 


eff 


1/2 


0.028 


(34) 


For the present airplane the wing stiffness in fundamental bending 
at 11.7 cps is approximately 19,000 pounds per foot. This value for k 
together with the dimensions given in table I and the estimate of equa- 
tion (34) gives the following relation for the root-mean-square buffeting 
shear at the root of each wing panel: 
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(35) 

and for the maximum buffeting shear likely to he encountered in a time At 

(eq. (B33)): 




/q iog e (n.7 At) 


62 


( 36 ) 


The least-squares relationship for the wing loads of the present 
tests with the basic airplane, equation (26), gives as a limit for very 
abrupt stalls, J 


/q log e (li.7 At) 


65.6 


while for very gradual stalls the limit is 


(37) 




= 3k 


/q loge(ll .7 At) 

and for the data as a whole, equation ( 4 ) and table V, an average is 


(38) 


/q. iog e (n.7 At) 


= 44.4 


(39) 


In view of the limited knowledge available about buffeting as a 
stationary random process, the number and character of the assumptions 
in appendix B, and the limited applicable experimental data on the aero- 
dynamic characteristics of stalled wings, the agreement between the 
constant value 62 of equation (36) and the values 65.6, 34, and 44.4 
obtained by least squares (eqs. (37)3 ( 38 ), and (39)) maybe fortuitous. 
The agreement shown does suggest, however, that further investigation 
is warranted of both the aerodynamic parameters and their relationship 
to the buffeting of other airplanes. 
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Buffeting coefficients .- The results of the present tests indicate 
that the usual buffeting coefficient of the 'form AL/qS would, for both 
wing and tail loads, be overly conservative if coefficients based on 
loads measurements at high altitudes were used for the estimation of 
loads at low altitudes. The tests also indicate that, for a given 
airplane, a simple comparison of l oads on the b asis of values of the 
dimensional forms AL/\/q or AI</(/q log e (f n At) would give more con- 
sistent results. To the extent that the simplified analysis of 
appendix B represents the buffeting of a straight-wing airplane in 


stalls, a coefficient of the form AT. 


//qkcs(l - e" A / 2 )log e (f n At) 
» A/2 


would be required to include both the geometry and the elastic properties 
of the wing, as well as the operating conditions of speed and altitude. 
Such a coefficient for the present abrupt- stall data would have a v al ue 
of approximately 0.03. Whether such a coefficient established for one 
type of airplane would give useful information about another type dif- 
fering, say, in wing thickness ratio or airfoil section would depend on 
the aerodynamic characteristics of the wing in stalls, as represented 

~/(Cla) € 


in the term c n £ 


Jeff' I 11 the absence of more experimental data 
on a spectrum of aerodynamic excitation for buffeting and on the effects 
of Mach number and angle of attack on both the spectrum and the aero- 
dynamic damping, a conclusion about a final form of a wing buffeting 
coefficient cannot be reached. However, should the results for the 
present unswept-wing airplane be confirmed for other similar airplane 
types, it should be possible to extend them to swept wings and to tails. 


Comparison of wing buffeting loads and design loads .- The results 
of the least-squares analysis of the wing buffeting loads of the present 
tests can be used to compare the maximum wing buffeting loads likely to 
be encountered in stalls with the wing design loads for the North 
American F-51D airplane. From equation (2 6 ) the amplitude of the maxi- 
mum buffeting- load increment in an abrupt stall of duration At is 
approximately 


Altf = 65 . 6 ^. log e (ll .7 At) 


(bO) 


The dynamic pressure of the stall can be expressed in terms of load 
factor, wing loading, and airplane normal-force coefficient as 

_ n(W/S) 
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Therefore AL^ can also "be expressed as 


Altf = 65.6 B&M log e (ll.7 At) 


(^1) 


The largest value of AI^ would be found in stalls at limit load factor 
at such speed and altitude that Cjfgg is as small as possible. The 

least value for Cu BB in stalls, figure 4(a), is 1.04. The limit load 

factor for the test airplane is 7.1 for a gross weight of 9,000 pounds. 
These values give, for the maximum value of AT W expected, 

^Wmax = 1^050^1og e (ll.7 At) (42) 

or, for a stall of 5 seconds' duration, ALu^ = 2,650 pounds. 

Such a buffeting load encountered in a stall at limit load factor 
would be superimposed on a steady wing-panel root structural shear of 
approximately 22,000 pounds. In terms of a gross weight of 9,000 pounds, 
a root-shear fluctuation of ± 2,650 pounds corresponds to a load- factor 
fluctuation of approximately ± 0 . 30 . 

Fatigue .- For fatigue studies, information is needed on the number 
of times a given value of load is exceeded in a given period. For a 
stationary random process, this information is provided by the mean- 
square load and the power spectrum of the load, as in equation (B26). 

The simple buffeting model considered in appendix B is a single-degree- 
of-freedom system which is very lightly damped. For such a system, the 
response to a random input has the character of a sine wave with a fre- 
quency roughly equal to the system natural frequency and an amplitude 
which fluctuates irregularly. The irregular a m plitude fluctuations 
are characterized by the probability distribution of equation (B3l) 
which gives the number of peaks per second which will exceed a given 
value. Since the total number of positive peaks per second corresponds 
to the natural frequency of the system f^ (with an equal number of 
minimums), equation (B3l) provides a simple basis for considering the 
fatigue aspects of buffeting. (See also ref. 7 .) Altho ugh based on a 
simplified model which ignores any contribution of higher vibration 
modes to the wing buffeting loads, equation (B3l) may well represent a 
satisfactory engineering approximation since modes of frequency higher 
than first bending ordinarily make but a small contribution to wing- 
root shear. 
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CONCLUDING REMARKS 


Wing and tail buffeting loads have been measured on a fighter-type 
airplane during I9A maneuvers. The half-amplitude of the largest fluc- 
tuation in structural shear was used as the measure of buffeting 
intensity in each maneuver. Correlation coefficients of 0.9 were found 
for loads on the left and right wings and the left and right horizontal • 
stabilizers. Least-squares methods have been used to illustrate certain 
trends in the data; in these studies the loads in the stall regime were 
found to follow a pattern which differed from that found in the shock 
regime. 

In the stall regime primary variables affecting the magnitude of 
the loads were speed and altitude as represented by the dynamic pressure, 
but the square root of the dynamic pressure was a better measure of the 
load than was the first power, a result which may be due to the action 
of aerodynamic damping. The loads measured in maneuvers of longer dura- 
tion were, on the average, larger than those measured in man euvers of 
short duration, a result which is in accord with considerations of sta- 
tionary random processes. As compared with abrupt pull-ups, load 
alleviation of about 50 percent was obtained by a gradual entry into 

the stall. 

/ 

In the shock regime, the primary variable affecting the magnitude 
of the loads was the extent of the penetration beyond the buffet boundary. 
The data do not provide a clear indication of a dependency of load on 
dynamic pressure, a result which may be in part attributable to the 
operating limitations of the airplane which restricted the range of 
the investigation in the shock regime; a more detailed investigation 
appears to be required. 

loads were also measured on a modification of the airplane incor- 
porating internal wing-tip weights which reduced the natural frequency 
of the wing in fundamental bending from 11 . 7 to 9.3 cps. Analysis of 
the measured loads indicated a reduction in wing loads of about 15 per- 
cent, and a similar percentage increase in the tail loads, as compared 
with the loads on the basic airplane. 

The loads on a simplified wing buffeting model have been examined 
on the assumption that buffeting is the linear response of an aerodyna- 
mically damped elastic system to an aerodynamic excitation which is a 
stationary random process. The results of the present tests for stalls 
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are sufficiently consistent with the results of the analytical study 
to suggest the examination of the buffeting of other airplanes on the 
same basis. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 11, 1954. 
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APPENDIX A 

SUMMARY OF STATISTICAL PROCEDURES 


A typical problem in linear regression involving a dependent vari- 
able w and, say, two independent variables x and y, which is solved 
by least- squares methods, is usually represented as finding the unknown 
coefficients a, b, and c in the equation 

w = ax + by + c 

given a set of N values of x and y assumed to be exact, and N 
corresponding measured values of w denoted by w'. For any set of 
values of a, b, and c, each measured value w'^ and the corre- 
sponding calculated value wq differ by the residual e-^ where 


c i = w 'i - Vi 

= v'i - ax-L - by-L - c 

The theory of least squares assumes that the "best" values of a, b, 
and c are those for which the sum of the squares of the resid- 
N 

uals €j_ 2 is a minimum, a condition which is fulfilled by the 

i=l 

values of a, b, and c in the so-called least-squares normal equa- 
tions which may be represented in matrix form as 


r 


N 

z* 

Ir 


V 




Z * 2 

z*y 

A 

a 

>= < 


£y 

z*y 



b 




where the summation 


denotes ^ . The resulting plane ax + by + c 
i-1 


passes through the point (w',x,y) determined by the mean values of w', 
and y . 
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The present paper is concerned with the application of least- 
squares methods to equations of the type where c = 0, and 

w = ax 


w = ax + by 

that is, problems where the least-squares line or plane is required to 
pass through the origin (w = x = y = 0). In this case for two inde- 
pendent variables, x and y, the values of a and b are given by 
the normal equations 



The solution may conveniently be written in terms of the inverse matrix 
which for second-order matrices is given by 



X* 2 X xy 


X 2 

1 

X* 2 Xy 2 -(l*y) 2 & y 




Accordingly 
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The sum of the squares of the residuals is given by 


I* 2 - IV * - ^ 



A measure of the spread in the measured values of w' is s w ' , the 
standard error of w' defined by 


s w ' = 



where w' is the arithmetic mean of the measured values . The 

N 

standard error of w' is usually most easily evaluated by the equation 



The standard error of the mean s w ' is proportional to s w ' and 
inversely proportional to the square root of the number of points, 
that is, 


s w ’ = 


w 


t 


✓n 


A measure of the ability of the regression equation to represent 
the data is given by the standard error of estimate of the equation, 
which for w = ax is 
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and for w = ax + by is 





The standard errors of' estimate of the constants a and b are related 
to the standard error of estimate of the equation and the terms on the 
principal diagonal of the inverse of the matrix of the coefficients of 
the normal equation by the relations 


s a ~ / c n s v 


s b = 1^22 s w 


The standard error of w' , that is, s w ', is a measure of the error 
involved in representing the N values of w' by their mean value w 7 ". 
An equation, say, w = ax, for which the. standard error of estimate s w 
is smaller than Sy' would ordinarily be considered an improvement over 
the mean-value representation, since it implies that specification of a 
value of x gives better information about the value of w' than does 
the mean value w' . The methods of the analysis of variance give a 
statistical estimate of whether the equation w = ax is improved by 
the addition of another variable y, to give w = ap* + bpy. For this 

particular question (see ref. 8) if ^"eq 2 and ]>~ep 2 represent the 

sum of the squares of the residuals of the one- and two-parameter equa- 
tions being compared and the ratio 


F 


In 2 - X ^ 2 



exceeds a certain critical value, then, on the basis of the evidence at 
hand, the chances are at least 100 to 1 that the improvement is real. 
The magnitude of the critical value of F depends upon the number of 
values N. For N = 25, 50, and 100, the values of F are 7. 97, 7«20, 
and 6.91, respectively. 
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Although linear dependency between two variables w and x is 
usually expressed by a relationship of the type w = ax + c when the 
measured values of x are considered exact, or in any event more nearly 
under experimental control than the measurements of w, there are 
instances when a more general measure of the linear dependency of two 
variables is desired. The coefficient of linear correlation r is 
such a measure which does not depend on the choice of w or x as 
independent variable or on the units of w and x. The value of r 
is usually calculated from the relation 



|^[nXw 2 “ (Zw) J |nX * 2 ~ (Z x) j 


but it can be shown that this value is equal to the square root of the 
product of the slopes a and a' in the two regression equations 

w = ax + c 

and 


x = a'w + c ' 


that is 



The values of r fall within the range -1=0=1, unit values indi- 
cating exact linear dependence and zero indicating complete independence 
of the two variables. A negative correlation coefficient indicates 
inverse dependency; that is, increasing values of one variable are 
associated with decreasing values of the other. 

For convenience in computation, all of the summations required in 
regression and correlation studies of the variables w and x may be 
obtained by expressing the N pairs of related measurements such as 
(w,x)j_ in the rectangular matrix 
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and premultiplying this matrix by its transpose ||m T ||, so that the 
following symmetrical square matrix results 



Similar considerations apply, of course, to the study of w, x, and y. 
More detailed treatment of the precision and interpretation of regression 
studies will he found in references 8 and 9» Numerical procedures are 
described in references 10 and 11. 
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APPENDIX B 

LOADS ON A SIMPLIFIED WING BUFFETING MODEL 


References 2 and 5 have illustrated the application of methods 
developed in the study of stationary random processes^ to the problem 
of the buffeting of an elastic structure such as a tail located in a 
turbulent airstream. A simple parallel treatment is possible which 
illustrates the form of the relationship between the airfoil motions 
and pertinent structural, geometric, and aerodynamic parameters for 
an elastically restrained airfoil subjected to the excitation of its 
own separated flow. 


The simplified model considered in the present section is a rigid 
airfoil of mass m, span b, mean chord c, and area S restrained by 
a spring of stiffness k to oscillate in vertical motion only. The 
vertical displacement z(t) from equilibrium can be expressed by the 
differential equation for a single-degree-of-freedom system: 


d^z ~ dz 2 F(t) 
— o + 2 7 <D n 77 + z = 
dt 2 dt m 


(Bl) 


where 7 is the ratio of the damping to critical damping, is the 

undamped natural circular frequency given by the relation 

(D n 2 = if (B2) 

and F (t ) is an impressed force. For an airfoil in a stream of air of 
dynamic pressure q, the exciting force associated with a time-varying 
fluctuating section normal- force coefficient c n (t) would be (three- 
dimensional effects being ignored) 

F(t) = c n (t)cbq (B3) 


^Time variations of a quantity during a particular time interval may 
be studied by the method of Fourier analysis, and this method can be gen- 
eralized to apply to a continuing nonperiodic disturbance through use of 
the concept of a stationary random process. This concept applies when 
the underlying physical mechanism which gives rise to an irregular dis- 
turbance does not change in time and the resultant process is thus both 
stationary and random. As a random process it can be described by cer- 
tain statistical parameters (mean, mean square, and power spectrum are 
ordinarily of chief interest); as a stationary random process, these 
parameters do not change in time and prediction on a statistical basis 
is therefore possible. For a more complete discussion see references 12 
and 13. 
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^ c n(t) is a random function of time but is expressible in terms of 
a power spectrum of the coefficient of the section-normal- force fluc- 
tuations c n (o>) such that the mean-square section normal-force coef- 
ficient is 



(Bl+) 


then z(t) is also a random function of time, expressible by a power 
spectrum z 2 (cd) and, by reason of equation (Bl), z 2 (co) is related 
to c n 2 (o>) through the admittance A 2 (cu) of the system by the relation 


C 2 b 2 n 2 

z2 (“) = - p - ' r c n 2 (ai)A 2 ((u) 
m ^°n . 

where 


(B5) 



The mean-square displacement of the airfoil is given by the definite 
integral of equation (B 5 ), that is. 


2 c 2 b 2 q 2 r 2/ \ 2/ V 
z = / c n 2 (o))A 2 (co) dn> (B 7 ) 

m^on 4 J 0 


Evaluation of the integral in equation (B7 ) could be a complex 
problem, evert under the assumption of positive damping, but, for small 
values of the damping, the admittance A 2 (oj) in equation (B 7 ) changes 
very rapidly in the frequency band in the vicinity of resonance, cu = u>n, 
and it is possible to substitute for c n 2 (u>) in equation (B 7 ) its 
value at and to write the approximate relation 


o c 2 b 2 q 2 p 
z^ ~ - c n ^ 


m' 


W 1 


(<%) / A 2 (o>) dca 


(B 8 ) 
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For the admittance given by equation (b 6), the area under the admittance 
curve is inversely proportional to the damping ratio since 



(B9) 


Therefore, the mean-square displacement is 


z 2 ^ Jtc 2 b 2 q 2 


Cn 2 (<Pn) 


(BIO) 


For the simplified buffeting model considered, aerodynamic damping 
forces would originate in the velocity of the vertical motion z and 
the damping ratio could be expressed as 


7 = 


qcb 

2nx%V 



(Bll) 


Where (pl^) e ff Be considered as an effective slope of the lift 

curve applicable to the damping of small bending motions of a stalled 
airfoil. The present flight tests have been concerned with values of 
wing root shear, which are analogous not to the airfoil displacement 
but to the load L = kz exerted on the spring support. Hence, an 
expression for the mean-square shear load in buffeting obtained from 

equations (B2), (BIO), and (Bll) would be L 2 = k 2 z 2 or 


rtkcbqV 


! ( c W) 


c n 2 (<%) 


eff 


(B12) 


Two characteristics pertinent to the definition of the spec- 
trum c n 2 (o>) are its level, as determined by the mean square, and its 
shape, or the frequency distribution of the excitation. These charac- 
teristics may be expressed by writing c n 2 (a>) in the form 


(B13) 


c n 2 (<J0) = c n 2 $(u>) 
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where $(<u) is the power-spectral-density function or shape parameter 
which defines the contribution to c n 2 from the excitation in any 
frequency band between o> and to + do>. Thus, in view of equation (b4). 



(B14) 


For a section property, it seems probable that the frequency cu 
is a less fundamental variable for defining the shape of the spectrum 
than a reduced frequency based on the speed V and a linear dimension 
related to the size of the airfoil or the chordwise extent of separa- 
tion. For the chord as the pertinent linear dimension, a reduced shape 

parameter is related to <t>(o)) by • requirements . of dimensional 

consistency, that is 



(B15) 


where the constant K which appears in the denominator is the area 
under the curve defined by the reduced shape parameter. Thus, on the 
basis of dimensional considerations, the spectrum c n ^(oo) may be 
written as 

=n 2 M - ^ i ®(f ) (B16) 

where 

k - r 0 kt) if) (bit) 


and the intensity of the fluctuations of section normal force at a 
particu lar frequency is seen to depend not only on the mean-square 

value Cn but also on the scale and speed and on the spectral dis- 
tribution of the excitation as expressed by the reduced shape parameter. 
From equation (Bl6), which provides a value for c n 2 (o>n), the mean- 
square buffeting load is 



jtkc^bq 


2K ( Cl u)eff 




(Bl8) 
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Little information is available concerning the shape parameter $ 

for stalled airfoils. In references 2 and 14 isotropic turbulence has 
been used to illustrate a random excitation expressible by a power 
spectrum. At a point in isotropic turbulence, the turbulent component 
of velocity w(t) normal to the free-stream velocity V results in 

+ ) 

an equivalent fluctuating angle of attack <x(t) = - which has a 

mean-square value and a spectrum a 2 (u>) that can be written in 

terms of a reduced frequency Zco/V as 



<x 2 (o)) 


= a 2 — 0 
KV 



(B19) 


where l is a linear dimension characteristic of the scale of the 
turbulence, and 


4 > 



3 Z 2 cd 2 



(B20) 


for which the constant K of equation (BI7) is equal to it. This 
particular shape parameter, which has been plotted in figure 19, is 
relatively constant and close to unity for values of reduced frequency 
less than 1 and then falls rapidly to low values. The assumption that 
the spectrum of the coefficient of the section-normal- force fluctuations 


( C0)\ 

—j has a shape similar to that expressed in 

equation (B20) with Z = c leads to an estimate of jt for the con- 
stant K in equation (Bl8) and provides a guide for estimating the 

(?)• 


on a stalled airfoil 0 


value of 01 


In equation (B3) and thus in equation (Bl8), section properties 
have been applied to the excitation of the entire wing, an application 
which, in general, would be expected to overestimate the net excitation 
since fluctuations at one chord station would not necessarily be in 
phase with fluctuations at another station. A simple overall correc- 
tion is possible, however, which is based on a correlation function 
observed in isotropic turbulence and is directly related to the spectrum, 



40 


NACA TN 3080 


equation (B19). This correction is similar to the length correction 
used in hot-wire anemometry and is used in reference 14 to relate the 
mean-square angle-of-attack fluctuation at a point along the span to 
the mean-square value over the entire span. It involves the ratio of 
the scale of the turbulence to the span b. If the same overall cor- 
rection is applied to the coefficient of section-normal- force fluctua- 
tions to t ake care of the major effects of spanwise load correlation, 

the wing C^ 2 would be related to the section c n 2 by the equation 

C N 2 = |(l - e“ b//c ) (B2l) 

This same overall correction leads to the final expression, applicable 

to the simplified model, for L 2 the mean-square force exerted on the 
model support 


.2 _ kc^bq 1 - e 


. e-Vo 




<S) 


■O' 


eff 


(B22) 


With slight modification, an expression applicable to the root 
shear of a wing panel can be obtained from equation (B22). For wing 
motions which are simplified in that only fundamental bending at natural 
frequency 0 % is considered, the vertical motion varies along the semi- 
span direction y in accordance with the shape of the bending mode zq(y) 
(taken as unity at the tip). The stiffness k would be an effective 
stiffness corresponding to this mode, where 

k = niecun 2 (B23) 

and me is an effective mass for bending in this mode, given by the 
integral of the product of the spanwise wing mass distribution m(y) 
and the square of the mode shape, or 

r b/2 ? 

1 % = / mCyjzp^y) dy 
J 0 


(B24) 
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Thus for the assumed wing, the mean-square root buffeting shear for one 
wing panel of span b/2 would be 


o kc^qb 1 - e 
L — — — 


-b/2c 


-n 


b /® ( c O 


eff 


V v / 


(B25) 


or, in terms of aspect ratio A = and total wing area S = be, the 

c 

me an- square root shear would be 


2 kcSq 1 - e 


-A/2 


L *> 


"n 


A/2 


( c O 


eff 


$ ( v/ 


(B26) 


For a given structure (c and fixed) the proportionality between 

I ?- and q (or V^) could be modified by changes in the value of the 

shape parameter with speed. If, however, the value of the 

reduced frequency lies in a nearly flat portion at the low- 

frequency end of the spectrum, then for a spectrum with a shape param- 
eter like that given by equation (B20), the value of the shape param- 
eter in equation (B26) can be replaced by its approximate value, 

unity, and 


, -A/2 n 2 

g «, MS j— 1 , 

4 ' ' 


A /2 K) 


(B27) 


eff 


Such a substitution would be valid over a range of speeds which is 
wider for low values of c and low values of natural frequency opn. 

The foregoing development deals with the mean-square load on a 
wing panel. If the buffeting of the simplified model can be considered 
a normally distributed stationary random process, then the relationship 
between the mean-square root shear L 2 and the probable amplitude AL 
of the maximum fluctuation occurring in a time interval At is fixed 
by the power spectrum of the load L^(cu). As shown in reference 12, 

(in the notation of the present paper) the number of peak values per 
second which will exceed a particular level ALp is approximately, 
when ALp is large, 
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poo 

/ o) 2 L 2 (<jo) dtu 

J n 


-+ 1/2 


j noo 

L 2 (co) dn> 
n 


-AL! 2 /2L 2 


(B28) 


Just as equation (B7) was simplified to equation (b 8) the term in 
brackets is easily evaluated, since 





and, for an admittance given by equation (b 6), 




d+jo 


(B29) 


(B30) 


Therefore, since = 2itf n , 


-AL 1 2 /2L 2 

n AL]_ = f n e 


(B31) 


and a value of AL will, on the average, be exceeded once in a time 
interval At given by the expression 


— = f 
At 


-AL 2 /2I? 


n c 


(B32) 


or the value AL which occurs once, on the average, in a time 
interval At is given by the equation 


AL 


i 


2L 2 logg ( f n At ) 


(B33) 
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The ratio AL/]/l 2 is plotted in figure 20 for two values of f n , 
9.3 and II.7 cps , corresponding to the basic and modified wing in the 
fundamental bending mode, the predominant mode in the wing buffeting 
time histories observed in the present investigation. 

Combination of equations (B27) and (B33) leads to an equation which 
relates the maximum load AI# (as measured in the present tests) in a 
stall of duration At to the geometric, structural, and aerodynamic 
characteristics of the simplified wing, 


at . kcS 1 
AIw * ( — 2 


-A/2 


1/2 r=- 


A/2 


n 


(%) 


\j<l loge(f n At) (B34) 


eff 
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TABLE I.- GEOMETRICAL DATA FOR TEST AIRPLANE 


Wing: 

Span, ft • . . 37.03 

Area, sq ft 240.1 

Mean aerodynamic chord, ft 6.63 

Aspect ratio 5*71 

Root thickness ratio 0 .I 5 

Tip thickness ratio 0.12 

Taper ratio 0.462 

Horizontal tail: 

Span, ft 13.18 

Area, sq ft 4l.O 

Weight at take-off, lb: 

Basic airplane 8,995 

Modified airplane 9,149 

Center-of-gravity position at take-off, percent M.A.C.: 

Basic airplane . . . 27.2 

Modified airplane 25.3 
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TABLE II.- NATURAL FREQUENCY OF AIRPLANE COMPONENTS 


Basic Modified 

airplane airplane 


Wing: 

Fundamental bending frequency, cps 11 • 7 9*3 

First asymmetric bending frequency, cps ... 22.3 18.1 

Torsion frequency, cps 38.0 34.5 

Second symmetric bending frequency, cps ... 5^.0 

Horizontal stabilizer: 

Primary bending frequency, cps 25 . 0 25 . 0 

First asymmetric bending frequency, cps ... 36 . 0 36.0 

Torsion frequency, cps 70*0 JO.O 


Fuselage: 

Torsion frequency, cps .... 
Side bending frequency, cps 
Vertical bending frequency, cps 


9.8 

12.5 

14.9 


9-8 

12.5 

14.9 



TABLE III.- BUFFETING CONDITIONS AND LOADS - BASIC AIRPLANE 
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Q LetterB used in this column have the following significance: 

vl very light buffeting L left roll-off 

I light buffeting R right roll-off 

m moderate buffeting N no roll 

h heavy buffeting 
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^ee footnote 



TABLE III.- BUFFETING CONDITIONS AND LOADS - BASIC AIRPLANE - Continued 
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TABLE III.- BUFFETING CONDITIONS AND LOADS - BASIC AIRPLANE - Continued 
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TABLE IV.- BUFFETING CONDITIONS AND LOADS - MODIFIED AIRPLANE 
(a) Operating Conditions 


Run 

mb 

c nb 

hp, 

ft 

n 

oc 

V 

At, 

sec 


c n e 

Pilot’s 

note 

(a) 

Pull-ups at an altitude of 30,000 ft 

1 

0.275 

1.194 

29,900 

0.10 

0.53 x 10" 5 

4.25 

0.236 

1.177 

Z, L, R 

2 

.287 

1.201 

29,900 

.50 

2.35 

1.90 

.273 

0 

vZ 

3 

.287 

1.197 

29,900 

.20 

.94 

8.12 

.249 

.984 

vl 

4 

.298 

1.250 

29,900 

.60 

2.51 

1.85 

.287 

.701 


5 

.327 

1.226 

29,800 

-50 

I.58 

3.80 

.298 

.601 

l 

6 

.354 

1.147 

30,300 

.60 

1.52 

. 2.90 

.327 

.765 

m 

7 

.377 

1.154 

29,900 

.80 

1.65 

3.50 

.324 

•791 

Z 

8 

.421 

1.083 

30,300 

1.10 

I.65 

3-53 

•357 

.728 

m 

9 

.452 

1.093 

30,400 

1.10 

1.36 

3.30 

.387 

.863 

m 

10 

.469 

1.094 

30,900 

1.10 

1.24 

4.13 

.368 

.815 

m 

11 

.474 

1.103 

29,800 

2.40 

2.48 

3.00 

.405 

• 524 

h, R, L 

12 

.476 

1.115 

30,200 

2.20 

2.28 

2.60 

.403 

.840 


13 

.479 

1.128 

30,100 

7.70 

7.80 

2.40 

.416 

•951 


14 

.480 

1.100 

30 , 4 oo 

1.60 

1.64 

2.60 

.414 

1.023 

m, R, L 

15 

.U83 

1.088 

30,200 

1.50 

1.49 

2.70 

.381 

.902 

h, R, L 

16 

.490 

1.092 

31,400 

2.20 

2.22 

4.00 

.399 

.848 

m 

17 

. 5 H 

1.107 

29,600 

2.80 

1.17 

2.20 

.469 

.783 

h 

18 

.513 

1.085 

29,200 

1.70 

.69 

3-20 

.443 

.644 

h 

19 

.516 

1.080 

30,300 

.90 

•75 

2.00 

.465 

1.050 

m 

20 

.516 

1.072 

30,300 

1.00 

.90 

2.70 

.466 

.758 

h 

21 

•537 

1.078 

30,700 

1.40 

1.04 

2.40 

.493 

.858 

h, R, L 

22 

.541 

1.072 

30,900 

1.00 

•73 

4.2 

.436 

.941 

d, R , L 

23 

.548 

1.087 

31,900 

2.10 

1.55 

2.00 

.505 

.648 

m, R, L 

24 

.716 

.629 

28,800 

— 


1.50 

.691 

.705 

Z 

25 

.743 

.512 

28,500 

— 


3.10 

.708 

.667 


26 

.758 

.400 

29,000 






.696 

.739 


27 

.763 

.332 

28,200 

— 


5.70 

•724 

.580 

m 

28 

.768 

.385 

29,000 

— 



.708 

.671 

h 

29 

.769 

.423 

29,300 

— 


5.40 

.697 

.720 

l 

30 

.771 

.340 

28,600 

— 


7.00 

.705 

.690 

m 

31 

.771 

.429 

29,700 




5.50 

.695 

.677 ! 

Z 

32 

• 773 

.279 

28,900 

— 


— 

.696 

.7^0 

l 

33 

•777 

• 334 

28,300 

— 


6.50 

•704 

.720 

Z, R 

34 

•796 

.273 

28,700 

— 


— 

.716 

.638 

Z 




Turn 

at an altitude of 30,000 ft 




33 

1 

0.350 

1.170 

28,000 ~ 

0.10 

0.23 X 10-3 

1.60 

0.349 

0.977 

R 




Pull-ups at an 

altitude of 10,000 ft 




36 

O.168 

1.241 

10,100 

0 

0 

4.00 

0.174 

0.699 

vZ 

37 

.205 

1.308 

10,200 

.40 

2.06 X 10-3 

3.40 

.178 

.904 ! 

Z 

38 

.207 

1.482 

9,950 

1.50 

7.48 

2.40 

.184 

.620 

Z 

39 

.261 

1.284 

10,100 

.70 

1.75 

1.60 

.232 

.862 

m 

40 

.267 

1.435 

9,950 

2.00 

4.63 

1.70 

.241 

.821 

h 

41 

.267 

1.593 

10,150 

4.00 

9.33 

1.55 

.247 

.468 

h 

42 

.332 

1.229 

10,700 

.70 

.87 

1.35 

• 299 

.781 

h 

43 

• 337 

1.271 

10,400 

1.60 

1.88 

2.22 

.280 

.700 


44 

.342 

1.338 

10,100 

7.60 

8.44 

1.15 

.329 

.298 

h 


a Letters used in this column have the following significance: 


vZ 

very light buffeting 

L 

left roll-off 

Z 

light buffeting 

R 

right roll-off 

m 

moderate buffeting 

N 

no' roll 

h 

heavy buffeting 
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TABLE IV.- BUFFETING CONDITIONS AND LOADS - MODIFIED AIRPLANE - Continued 
(b) Left and Right Wing Loads 



Left 

Right 


M 

Ib/sq ft 

Cn 

ac n 

At load/ 
sec 

AL, 

lb 

M 

q, 

lb/ sq ft 

cn 

ac n 



Pull 

-ups at 

an altitui 

ie of 3 

0,000 f 

t 





1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 
22 

23 

24 
23 

26 

27 

28 

29 

30 

31 

32 

33 

34 


445 

405 

300 

445 

445 

386 

500 

560 

613 

695 

605 

640 

810 

735 

730 

600 

965 

800 

710 

650 

607 

825 

715 

272 

355 

325 

214 

220 

440 

382 

805 

795 

228 

600 


0.261 

.271 

.277 

.287 

.323 

.328 

.351 

.384 

.402 

.*440 

.425 

.453 

.458 

.468 

.467 

.454 

.487 

.450 

.497 

.493 

• 530 

.501 

.516 

.714 

.718 

.752 

.724 

.721 

.721 

.763 

.7^7 

.770 

.7V7 

• 753 


30 

32 

33 
36 
46 

46 

54 

64 

68 

81 

79 

89 

91 
94 
94 

84 

106 

92 
108 
105 

110 

105 

106 

236 
243 

267 

260 

256 

237 
280 

252 

283 

269 

278 


1.208 

1.050 

1.007 

.857 

1.175 

.890 

1.146 

.861 

.920 

.894 

1.025 

1.067 

.983 

1.048 

1.093 

1.001 

1.020 

1.109 

.967 

.957 

1.010 

.966 

.971 

.651 

.657 

.513 

•579 

.607 

.660 

.451 

.638 

.478 

.588 

.583 


.030 

.080 

-.012 

-.002 

.051 

.076 

.175 

.148 

.125 

.156 


1.98 

1.25 

1.92 

1.54 

1.92 

2.45 

1.85 

2.09 
2.8o 
1.40 

1.70 

.76 

• 75 
.56 
.54 

1.42 

1.18 

2.59 
.70 

1.28 

1.50 

1.76 

1.39 

.20 

2.48 

2.44 

5.70 

6.57 

4.01 

3.10 

2.32 

2.20 

5.60 
5.72 


420 
500 
67O 
370 

344 

575 

540 

421 

705 

775 

785 

745 

610 

700 

640 

820 

725 

845 

560 

543 

600 

595 

140 

245 

455 

265 

295 

430 

142 

620 

650 

227 

510 


.278 

.282 

.294 

.318 

• 349 
.357 

• 394 
.432 
.441 

.417 

.458 

.463 

.448 

.431 

.434 

.481 

.449 

•497 

.492 

•550 

.502 

.520 

•714 

•730 

.756 

.7^0 

.752 

-759 

.763 

.7^6 

.75^ 

-759 

.767 


30 

1.198 



34 

1.341 


35 

1.190 


58 

1.248 


44 

1.090 


53 

1.140 1 


56 

1.073 


67 

1.056 


80 

1.112 


81 

1.022 


76 

•973 


90 

1.068 


93 

1.005 


86 

1.065 


79 

1.102 


76 

.988 


104 

1.106 


91- 

1.06l 


108 

.966 


104 

.949 


112 

1.001 


105 

.960 


108 

.977 


236 

.662 

0.011 

251 

.584 

.029 

267 

.487 

.075 

272 

.586 

.086 

279 

.421 

-.015 

263 

.540 

-145 

280 

.452 

.080 

251 

.635 

.166 

275 

.578 

.157 

278 

.513 

.118 

287 

.520 

.169 


1. 

.32 

.60 

-38 

1.20 

.60 

1.42 

1.54 
.80 

1.34 

2.02 

•59 

.60 

1.20 

I.65 

2.13 
1.45 
2.71 

-70 

1.34 

1.30 

1.75 

1.24 

.20 

1.64 

1.14 
4.70 

3.32 

1.08 

3.10 

2.39 

4.54 

2.60 

4.61 






Tu 

m at an 

altitude 

of 30, 

000 ft 



35 

48o 

0.350 

59 

1.156 

1 

1 

1 

1 

1 

1 

0.4o 

391 

0.350 

59 

1.145 


o.4o 



Pull 

-ups at 

an altitut 

Ie of 1 

p 

Hi 

t 





36 

37 

38 

39 

40 

41 

42 

43 

44 


315 

510 

785 

655 

710 

850 

450 

905 

1,380 


0.165 

.196 

.200 

.246 

-253 

.257 

.319 

.327 

.336 


27 

38 

4i 

61 

65 

67 

100 

90 

114 


1.273 

1.254 

1.411 

1.150 

1.082 

1.130 

1.225 

1.028 

1.150 


2.22 

1.01 

.62 

.80 

• 91 

.69 

.50 

1.39 

• 39 


435 

585 

510 

810 

830 

805 

1,450 

1,085 

1,210 


O.165 

.187 

.199 

.247 

.261 

.257 

.313 

.316 

.337 


27 

36 

40 

61 

69 

67 

96 

99 

114 


1.180 

1.194 

1.174 
1.150 

1.175 

1.130 

1.182 

1.120 

1.160 


2.76 

1.88 

.69 

•79 

.46 

.72 

.76 

.83 

.38 
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TABLE IV.- BUFFETING CONDITIONS AND LOADS - MODIFIED AIRPLANE - Concluded 
(c) Left and Right Tail Loads 



36 

37 

38 

39 

40 

41 

42 

43 

44 


295 0.165 27 1.290 2.04 

435 .190 36 1.218 1.65 

505 .200 41 1.410 .62 

540 . 247 61 1.160 .74 

520 .262 67. 1.158 .61 

595 - 263 70 1.380 ..30 

495 <309 86 1.100 .94 

700 .324 104 1.198 .52 

915 .334 113 .910 .64 


310 0.165 27 1.180 2.77 

517 .194 37 .898 1.20 

660 .199 40 I.380 .69 

440 .253 65 1.210 .43 

64 0 .259 68 1.147. .56 

715 -262 69 1.320 .39 

430 .311 95 1.160 .82 

595 -320 101 1.160 .69 

715 .336 113 1.080 .45 
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TABLE V.- SUMMARY OF WING- LOAD ANALYSIS - STALL REGIME 


Equation 

number 

Equation 

Sum of 
squares of 
residuals, 

I‘ 2 

Standard 
error of 
estimate, 
s, lb 

(1) 

AL^ = a^_ 

585 X 10^ 

255 

(2) 

AL^ = agq 

770 

293 

(5) 

= a 5 ^/q 

46l 

226 

00 

Alty = a^/q logg(ll.7 At) 

386 

206 

(5) 

f , -aZ/O.OCkV) 

AL^ = ya^ + b^e /y q 

341 

196 

(6) 

% = (a 6 + bge"^/ 0 * 0047 )^ log e (ll.7 At) 

287 

178 


Constants 


ai = 749 ± 27 
a 2 = 6.54 ± 0.27 
aj = 74.4 ± 2.4 
ajj. = 44.4 ± 1.3 
a^ = 111.5 ± 6.9' 
ag = 65.6 ± 3.8 


b 5 = -55-1 t 9-9 
bg = -31.6 ± 5.4 









58 


NACA TN 5080 


TABLE VI.- SUMMARY OF TAIL- LOAD ANALYSIS - STALL REGIME 


Equation 

number 

Equation 

Sum of 
squares of 
residuals , 

Standard 
error of 
estimate , 

S y lb 

(7) 

AI/j — A 

304 x 10^ 

184 

(8) 

ALj = Agq 

384 

207 

(9) 

1 

II 

280 

176 

(10) 

ALr = A 10 A log e (ll.7 At ) 

257 

.170 

(11) 

ALj = (An + b^^ 0,0 * 7 )^ 

174. 

140 

(12) 

ALp = (A12 + B 12 e < ^ C / 0 * 0 °4 V ^ (1 i oge (n.7 At) 

161 

135 


Constants . 


Af = 4l4 ± 19 
A Q = 3.59 ± 0.19 
A9 = 41.0 ± 1.8 
A]_ 0 = 24.4 * 1.0 
An = 75-4 ± 3-5 
A i2 = 44.1 * 2.9 


Bn = -51.2 ± 5.0 
b 12 = -29.2 ± 4.1 
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TABLE VII.- SUMMARY OF WING- LOAD ANALYSIS - SHOCK REGIME 


Equation 

number 

Equation 

Sum of squares 
of residuals, 

Standard error 
of estimate, 
s, lb 

(14) 

ALw = a^ 

834 X 10^ 

578 

(15) 

Alty = a^q 

1,224 

715 

(16) 

Aly = a l6 /q 

| 

1,009 

648 

(IT) 

AL w = a l7 P 

133 

238 

(18) 

ALw = a 1 0Pq 

192 

283 

(19) 

ALtf = a^P^q 

125 

228 


Constants 


a lA = 94 0 ± ll6 
a 15 = 2.8l + 0.28 
a x 6 = 52.2 ± 7.3 
a yi = 2500 ± 107 
a]_8 = 9.68 ± O.51 
a 19 = 155-5 ± 6.4 
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TABLE VIII.- SUMMARY OF TAIL- LOAD ANALYSIS - SHOCK REGIME 


Equation 

number 

Equation 

Sum of squares 
of residuals, 

Standard error , 
of estimate, 
s, lb 

(20) 

Alq. = A 20 

218 x 10 ^ 

295 

(21) 

ALp = A2]_q. 

334 

365 

(22) 

ALp = A 22 ^5 

270 

535 

(23) 

ALp = Ag^P 

67 

167 

(24) 

ALp = A24 p q. 

71 

173 

(25) 

ALp = 

73 

174 


Constants 


a 20 = 508 ± 59 
a 21 = 1*52 ± 0.24 
A 2 2 = 28.2 + 3.8 

A23 = 1254 ± 76 
A 2 4 = 4.59 ± 0.28 
A25 = 75.2 ± 4.6 



Posit: 

wing-1 

weighi 


Figure ] 
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Figure 2.- Side view of test airplane. 
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(a) Altitude, 30,000 feet. 



(b) Altitude, 25,000 feet. 

Figure 7.- Wing buffeting loads - basic airplane 
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Figure 11.- Correlation between left and right wing buffeting loads - 

basic airplane. 





Figure 13.- Correlation "between wing and tail buffeting loads - basic 

airplane . 
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Figure 15. - Relationship between buffeting load, ^ and Mach number in shock regime 
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Figure 1 6 .- Comparison of results of least-squares analysis with trends 
shown by method of averages . Circles represent data from figure 14 . 
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